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We investigate the oxygen octahedral rotations that occur in two perovskites, SrTiOa and PbTiOs, 
as a function of applied three-dimensional electric displacement field, allowing us to map out the 
phase diagram of rotations in both the paraelectric and ferroelectric regions of the polar response. 
First-principles calculations at fixed electric displacement field are used to extract parameters of 
a Landau-Devonshire model that is analyzed to identify the phase boundaries between different 
rotational states. The calculations reveal a rich phase diagram of rotations versus applied field in 
both SrTiOs and PbTiOa, although the details are quite different in the two cases. 
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I. INTRODUCTION 

Oxygen octahedral rotations can have a significant 
impact on the behavior of ABO3 perovskites, affecting 
electronic, dielectric, ferroelectric and magnetic prop- 
erties. For example, the octahedral rotations cou- 
ple strongly with the magnetic structure in transition- 
metal perovskites by modifying the metal-oxygen-metal 
bond angles that are critical to determine the magnetic 
interactions In some materials, such as SrTiOa, the 
oxygen rotations give rise to a non-polar antiferrodis- 
tortive (AFD) ground state contribute to the suppres- 
sion of ferroelectric (FE) order. However, the recent 
discovery of rotation-driven improper ferroelectricity in a 
superlattice^ has inspired a search for this type of ferro- 
electricity in other types of materials jii^r— Because they 
can also couple with magnetic properties, octahedra ro- 
tations offer a promising approach to the discovery or 
design of new multiferroic perovskites .iiSrJ^ 

In the cubic structure, SrTiOs and PbTiOa both show 
an AFD instability at the zone corner (i? point) of the 
Brillouin zoneJ^ii^ Following this instability leads to the 
experimentally observed tetragonal ground state with ro- 
tations along the [001] axis for SrTiOs, while PbTiOs 
prefers rotations along the [111] axis. For PbTiOs, how- 
ever, this is not the equilibrium structure; instead, a 
strong FE instability at T out-competes the AFD insta- 
bility, giving rise to a tetragonal FE ground state with- 
out rotations. Nevertheless, the AFD modes provide a 
potential source of instability in PbTiOs, as has been 
predicted for example for surfacei^iii and interface*' ge- 
ometries. (In SrTiOs, a weak FE instability at T is found 
in some calculations, depending sensitively on lattice con- 
stant and other details. However, experimentally the ma- 
terial just barely avoids this instability, remaining para- 
electric down to zero temperature.) 

While oxygen octahedral rotations clearly play an im- 
portant role in these and other perovskites, and while 
they are known to be strongly affected by stres g-'^^d^dQ 
and temperature^^'' we are not aware of any previous 
study of the effects of electric fields on the AFD rota- 
tions. In this paper, we study the phase transition be- 
havior of the AFD modes in SrTiOs and PbTiOs under 



three-dimensional constant electric displacement field. 

Of course, SrTiOs and PbTiOs are qualitatively dif- 
ferent in that the latter is ferroelectric while the former 
remains paraelectric down to zero temperature. However, 
the choice of fixed electric displacement field D for the 
boundary conditions in this study allows us to treat both 
materials on an equal footing. The situation would have 
been much more complicated if we had chosen to work 
at fixed electric field £, because the energy landscape is 
multivalued and the paraelectric configuration is unsta- 
ble at small 5 in a ferroelectric material like PbTiOs. At 
fixed D, however, the energy landscape remains single- 
valued, thus allowing access to the entire electric equation 
of state for PbTiOs as well as SrTiOs.^^'^^ Indeed, both 
materials have a large static dielectric constant, so that 
mapping at fixed D is qualitatively similar to mapping at 
fixed P. Thus, in our study, the main qualitative differ- 
ence between SrTiOs and PbTiOs will be related to the 
fact that the rotational instability that prefers to develop 
along a [001] axis in SrTiOs instead prefers a [111] axis 
in PbTiOs. 

The paper is organized as follows. In Sec. |IT]we intro- 
duce the Landau-Devonshire model, provide the details 
of our first-principles calculations, and specify the termi- 
nology for symmetries that will be used later. In Sec. IIIII 
we present the results of the first-principles calculations 
in one-dimensional D-field space and discuss the fitting 
of the model, which is then used to compute the phase di- 
agram of rotational phases in three-dimensional D space. 
Finally, Sec. II VI contains a summary. 

II. PRELIMINARIES 

A. Landau-Devonshire model 

In order to explore the octahedral-rotation phase dia- 
gram in the space of D fields, the internal energy U has 
to be calculated and minimized on a three-dimensional 
mesh of D values. Near the phase boundaries between 
different rotational phases, this process would be quite 
tedious; the first-principles calculations would need to 
be very carefully converged, and the procedure would 
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become quite time-consuming. We therefore introduce a 
Landau-Devonshire model to study the phase transitions 
in this system, with the coefficients in the model being 
obtained from fitting to our first-principles results on a 
smaller database of D values. This model can then be 
used to locate the phase boundaries efficiently. 

As mentioned above, the dominant AFD rotational in- 
stabilities for paraelectric SrTiOs and PbTiOa are both 
at the R point in the Brillouin zone (corresponding, in 
the most general case, to the a~b~c~ Glazer notation). 
We therefore focus on these modes here, and define a 



Ue{e) 

C/int(D,0) 

Here we have made the approximation of truncating the 
expansion systematically at overall fourth order, and 
a, /?, 7, /i, w, a, T, A and k are coefficients that need to be 
fitted from the first-principles calculations. The terms in 
(7int describe the coupling of D and octahedral rotations. 
There is no strain in this expansion since each term is 
defined assuming that the strain is fully relaxed for each 
(D, e) value. 

In the present work we are really only interested in the 
internal energy difference C/tot(D, 9) — Utot(D, 0) between 
the states with and without octahedral rotations. We 
denote this quantity simply as U and note that 

f/(D,0)==C/e(0) + C/int(D,0). (5) 

In order to fit the coefficients from first-principles cal- 
culations, we first apply D along just one Cartesian di- 
rection, which we choose as D^, to find the coefficients 
^,uj,a,T and A. For this case we set = Dy = and 
find 

u = ^le^ + TDlel + \Dl{el + el) 

+Loe^ + {a- 2Lo){elel + elel + elel). (6) 

We then do a series of calculations in which we choose dif- 
ferent initial structures with equillibrium rotations along 
0=[1OO], [110], [001] or [111] at I?^=O.Oa.u., and for each 
choice (and for each D^) we relax all the coordinates to 
obtain the internal energy U{Dz, 9). (We increase in 
increments of 0.04 a. u. up to 0.12 a. u. for SrTiOs, and in- 
crements of 0.02a.u. up to 0.08 a.u. for PbTiOa.) Fitting 
the model parameters to this first-principles database of 
information, we obtain all the coefficients in Eqs. ([3]|3]) 
except for k. We then do one more series of calcula- 
tions with both D and along the [111] direction, i.e., 



vector octahedral rotation = {0x,dy,dz) describing a 
rotation by angle 9x around the a;-axis, etc. (or, more 
generally, by angle = |0| about axis 9). 

Within our Landau-Devonshire model, then, the inter- 
nal energy Utot(D,d) is expanded as a function of dis- 
placement field D = {Dx, Dy, Dz) and octahedral rota- 
tions ~ {Ox, dy,9z) as 

C/tot = UD + Ug + C/i„t (1) 

where 



(2) 
(3) 

(4) 

I 

D = {Do, Do, Do) and = {0o,Oo,9o), for which the 
model predicts 

U = S/ifo + 3{oJ + (J)e^ + 3{t + k + 2X)T>lel. (7) 

(Here D^ is increased in steps of 0.02 a.u. up to 0.08 a.u. 
for both materials.) Fitting in a similar way to these re- 
sults, we obtain the parameter k as well. Once all the 
parameters are in hand, we can go back to Eq. ([5]) and 
study the full behavior of octahedral rotations as a func- 
tion of three-dimensional D space using this model. 

B. First- principles methodology 

Our calculations were performed within density- 
functional theory in the local-density approximation^^ 
using norm-conserving pseudopotentials^^ and a plane- 
wave cutoff of 60 Ha. A 6x6x6 Monkhorst-Pack grid32 
was used to sample the Brillouin zone. The unit cell for 
simulating the i?-point rotation is doubled to obtain a 
10-atom fee cell. The atomic coordinates and lattice vec- 
tors of this cell were relaxed until all atomic force com- 
ponents were smaller than 10~^ Ha/Bohr and all stress 
components were below 10^^ Ha/Bohr'^. We used the 
open-source ABINIT code packagers, with the implemen- 
tation of the constant-displacement-field method in 3- 
dimensions^^ to calculate the internal energy at a each 
specified D field. 

C. Terminology for symmetries 

Here we introduce the notations that we will use for 
describing rotational phases, following a similar scheme 



= a{Dl +Dl + Dl) + fi{Dt + Dl + D^) + 7(^'A' + DIdI + dIdI), 

= ^i{el + el + el) + u{et + e^ + el) + cj{elel + elel + opl), 
= T{elDl + elDl + olDl) + \{elDl + oIdI + elDl + oIdI + oIdI + elDl) 

+ K{9xeyDxDy + 9xezDxDz + e yO z D y D z) . 
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FIG. 1: Internal energy U of Eq. ([Sjl for D applied along the z direction, for phases with the octahedral rotations constrained 
to be about different axes as indicated in the legend, (a) SrTiOa; (b) PbTiOa. Symbols are from first-principles calculations; 
curves are from the Landau-Devonshire model. 



as the one often uses for polarization.^'^ When the octahe- 
dral rotation axis is constrained to a symmetry axis lying 
along (001), (111), or (Oil), the resulting phase becomes 
tetragonal (T), rhombohedral (7^), or orthorhombic (O), 
respectively Similarly, the M. phases arise when rota- 
tion axis is confined to a mirror plane. There are three 
cases: A^c, in which the axis is along [0,m, u]; and TWa 
or A^B, in which the axis is along [uuv] with u < v or 
u > V, respectively. The triclinic phase (Tri) occurs if the 
axis is along [uww] with u ^ v =/= w ^ 0. We also intro- 
duce the Cartesian subscript a = {x,y,z} to specify the 
unique Cartesian direction when needed. For example, 
Ta denotes the tetragonal phase with rotation axis along 
direction a, while Oa and A^ca denote the orthorhom- 
bic phase and A4c phases with rotation axis lying in the 
plane perpendicular to the a direction. Similarly, Mau 
and A^bq are the TWa and Mb phases with the non- 
equal component v in [uuv] along the a direction. 



III. RESULTS 

A. First-principles calculations 

We first carry out a series of calculations, starting from 
D = and increasing in steps of 0.04 a.u. for SrTiOa 
and 0.02 a.u. for PbTiOs, to explore the resulting be- 
havior for the case that the octahedral rotation is con- 
strained to lie along the [100], [110], or [001] axis. At 
each Dz, the structure is fully relaxed with respect to 
both ionic positions and lattice parameters. In all three 
cases in both materials, the rotations, which are fully 



developed at Dz = 0, are gradually suppressed with in- 
creasing Dz until they disappear completely at a critical 
value of Dz- We also attempt this procedure for the case 
that the octahedral rotations started along the [111] di- 
rection at D = 0. However, for SrTiOs a rotation along 
[111] is a saddle point, rather than a local minimum, of 
the D = energy landscape, and the breaking of the 
three- fold symmetry about [111] by the applied Dz im- 
mediately causes the rotation axis to switch to either the 
[110] or [001] direction. For PbTiOs, by contrast, the 
D = system has its minimum- energy AFD axis along 
[111], and we can also follow the evolution of this fourth 
case as Dz is applied. In this case we find that 9z grad- 
ually increases, and 9x = 0y gradually decrease, with 
increasing Dz, until a critical value is reached at which 
6x and By vanish and the solution merges with the one 
with the rotation axis constrained to [001]. The results 
of these calculations are shown as the symbols in Figs. [T] 
and m where the internal energy difference U of Eq. ([5]) 
and the equilibrium rotation angles are plotted versus 
Dz- 

We also carry out calculations for both materials with 
D and the rotation axis both constrained to lie along 
[111]. As mentioned in Sec. Ill Al the purpose of this is 
just to obtain the additional coefficient k that was not 
determined from the calculations with D along [001], so 
it was not necessary to study other rotation axes for 
this case. The results (not shown) again indicate that 
the rotations gradually decrease with increase of Dq for 
D — {Dq, Dq, Dq), although in the case of SrTiOs the 
rotations never vanish over the range of Dq up to 0.10 
a.u. studied here. 
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FIG. 2: Octahedral rotation angles {9x,9y,9z) for different 
phases for D applied along z. (a) SrTiOa; (b) PbTiOa. For 
[110] cases {9x=0y), 9^ is plotted. For the PbTiOs [111]- 
derived case {9x=9y^9z), dx and 9z are plotted. Symbols are 
from first-principles calculations; curves are from model. 



B. Fitting of the model parameters 



We now use the results of the above first-principles 
calculations to determine the parameters in Eqs. ([3]|4]) 
following the procedure detailed at the end of Sec. Ill Al 
The resulting parameter values are reported in Table HI 
The predictions of the fit (solid curves) are compared 
with the direct first-principles results (symbols) in Figs.[T] 
and [21 It is clear that the model agrees quite well with 
the first-principles calculations. 



TABLE 1: Fitted coefficients of the Landau-Devonshire model 
of Eqs. Pl5| . defined with energies in meV, rotation angles in 
degrees, and displacement fields in a.u. 









a 


r 


A 


K 


SrTiOa 
PbTiOa 


-0.863 
-0.661 


0.015 
0.019 


0.036 
0.033 


64.45 
78.80 


41.19 
96.27 


-139.36 
-147.91 



C. Details for D along [001] 

From Fig.lll we can see that SrTiOs and PbTiOs have 
different octahedral rotation patterns. At D = 0, SrTiOs 
has the lowest energy in the T phase, which is its true 
ground state experimentally below 105 K, and the highest 
energy in the TZ phase (not shown in the figure because it 
is destabilized by any finite Dz.) In PbTiOs, on the other 
hand, the energy ordering is just the opposite, with the 
TZ phase lowest and the T phase highest in energy. In the 
context of Eq.lHl the energy ordering of the phases at D = 
is determined by the combination of parameters (cr — 
2a;), with the TZ ot T phase lowest in energy when this 
combination is negative or positive, respectively. This is 
confirmed by the coefficients in Table HI 

As Dz increases. Fig. [T] shows that the internal energy 
U of SrTiOa and PbTiOs increases and finally reaches 
zero. Recalling that U is measured relative to the struc- 
ture with no rotations, we conclude that the octahedral 
rotations disappear at a sufficiently high Dz field. How- 
ever, the behavior is different for these two materials. 
Fig. [TJa) shows that for SrTiOs the T phase^^ with rota- 
tion axis along [100] or [010] has the lowest energy as Dz 
increases, while the T phase with rotation along [001] in- 
creases sharply in energy and becomes the least favorable 
state when Dz > 0.05 a.u. This suggests that the T phase 
with its rotation axis perpendicular to D is favored, while 
the one with rotation axis parallel to D is suppressed. 
However, PbTiOs has a quite different behavior, as can 
be seen in Fig. \V[h) . The state of lowest internal energy 
at D = is the TZ phase. As Dz increases, the rotational 
axis is perturbed to be along [uuv] for v > u, putting the 
system in the A4a phase. Eventually, the internal energy 
of this state merges into the curve for the T phase (axis 
along [001]), indicating a phase transition from A4a to T 
at some critical value of Dz- The T phase, with its rota- 
tion axis parallel to D, is then favored at higher Dz, until 
there is a second phase transition at which the rotations 
disappear. 

The details of the rotational behavior in PbTiOs can 
be seen more clearly in Fig. Hfb), which shows the vari- 
ation of the various rotation angles with Dz field. The 
rotation angles decrease as Dz increases for all phases ex- 
cept for the initial TZ phase. This phase is immediately 
perturbed to become A^a as soon as a non-zero Dz is 
present. With increasing Dz, the rotation angles 9^ and 
6y decrease, but 9z increases. That is, the rotation axis 
starts from [111] (TZ) and then rotates in the (lIO) plane 
(A^a) towards [001] (T). We can now see that the crit- 
ical Dz at which T phase is reached (i.e., at which 9x 
and 9y vanish) is at Dz=0-058a,.u. This also corresponds 
to the merger of A^a and T phases in the internal en- 
ergy curves of Fig. [U^b) . For larger Dz , 9z then decreases 
monotonically and reaches zero at £'2=0.092 a.u. 

For SrTiOs, on the other hand, the picture is sim- 
pler. As Dz increases in Fig. [^Ja), the rotation axis re- 
mains along [100] while the amplitude of 9^ monotoni- 
cally decreases and disappears entirely at a critical value 
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FIG. 3: Phase diagram for rotational phases of SrTiOs under applied D field. Each panel is a cut plotted in the Dx-Dy plane 
(note the change of scale between top and bottom panels) at fixed D^. (a) Dz— 0.00 a.u.; (b) 0.015 a.u.; (c) 0.08 a.u.; (d) 0.16 
a.u.; (e) 0.172 a.u.; (f) 0.24 a.u. Solid lines are second-order boundaries; decorated solid lines are first-order boundaries; and 
dashed, dash-dotted, and dotted lines are special cases of higher symmetry induced by high-symmetry D-vectors (see text for 
details). Gray area in (b) is detailed in Fig. ID 



of D^=0.144a.u. 



D. Three-dimensional D field 

We now turn to a detailed discussion of the behavior of 
SrTiOa and PbTiOs as a function of three-dimensional 
D field, based on the model of Eqs. using the co- 

efficients fitted from first principles as reported in Table 

m 

First, note that because the coefficient /i is negative 
in both SrTiOa and PbTiOs, we are guaranteed to get a 
phase with non-zero rotations at small D. Also, because 
of the non-zero value of k, we generically obtain a tri- 
clinic rotational axis {d^ Oy ^ d^) at a general point 
Dx ^ Dy ^ Dz ^ in D space. High-symmetry phases 
will only exist under special conditions, i.e., when one 
or more D components vanish, or when two or more D 
components are equal. 



1. Phase diagram for SrTiOz 

Figure [3] shows several two-dimensional {D^, Dy) slices 
of the three-dimensional phase diagram of STO taken at 



different values of Dz- In these panels, the outer solid 
(blue) boundaries (and also the inner ones in Fig. [3]^d- 
e)) indicate a second-order phase transition from a phase 
with octahedral rotations to a phase without rotations. 
Other solid lines represent first-order phase boundaries 
as will be explained below. Dashed and dotted lines 
are not true phase boundaries, but instead denote high- 
symmetry structures that occur as special cases along 
special lines or planes in D space; we use dashed lines 
(red) for A4a or Mb phases, dotted lines (black) for T 
phases, dashed-dotted lines (green) for Mc phases, and 
short-dashed-dotted lines (brown) for O phases, using the 
notation developed in Sec. Ill CI 

In Fig. [3lja), for — 0, the squarish sohd curve 
marked by crosses (green) is a first-order boundary sepa- 
rating the Tz phase (inside) from phases with Oz=0 (out- 
side). For generic {Dx,Dy) outside, this corresponds to 
the A4cz phase (recall, from Sec. Ill CI that this means 
that the rotation axis lies in the 6^ — 0y plane). Along 
the horizontal axis {Dy — 0) outside, shown by the dotted 
(black) line, the Tz and Ty phases are degenerate. How- 
ever, any small finite Dy favors the Ty phase and adds a 
small 6x component (via the k term) so that the Adcz 
phase results. That is, crossing this dotted line from 
negative to positive Dy just causes 9x to cross smoothly 
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FIG. 4: Enlargement of the central portion of Fig. [3] (Da:, Dy £ 
[-0.08, 0.08] a.u.) for small values, (a) D^=0.00 a.u.; (b) 
0.01 a.u.; (c) 0.015 a.u.; (d) 0.025 a.u.; (e) 0.04 a.u.; (f) 0.05 
a.u. Color coding is such that pure red, green and blue cor- 
respond to the Tx, Ty and Tz phases respectively, with color 
values weighted according to \8x\, \0y\ and \dz \ for intermedi- 
ate phases. 



through zero, so that this is not a true phase boundary. 

Along the [110] direction in D space, the behavior is 
rather complex. We let = Dy = D^. Recall that in- 
side the square region (small Dq) one finds the Tz phase. 
Next comes a segment of first-order phase boundary, in- 
dicated again by a solid line vifith crosses (green), along 
which there are two degenerate A4cz phases with rota- 
tion angles {9a,0h,0) and {9b,0afi) (with 6a ^ Oh). Any 
small step away from this line (while remaining in the 
Dz=0 plane) favors one or the other of these phases, and 
also slightly perturbs its angles 9^ and 9y. Thus, when 
crossing this line, both 6^ and 6y jump discontinuously. 
When Dq increases further, as shown by the dash-dotted 
(brown) line, one finds the Oz phase exactly on this line, 
but it is just a special case of the Mcz phase as \6x\ and 
\9y\ cross smoothly through each other. Like the dotted 
(black) line, therefore, this is not a true phase boundary. 

As Dz increases from zero, the behavior of the phase 



diagram is initially very complex, especially in the vicin- 
ity of the squarish central region of Fig. [2Jb). The phase 
behavior in the outer region at Dz — 0.015 a.u. is shown 
in Fig. ^h). The Mcz phase at Dz=0 is perturbed to 
become triclinic as 9z becomes non-zero linearly in Dz- 
The Tx and Ty lines at Dz—0 are perturbed to A^c struc- 
tures as shown by the dash-dotted (green) lines, and the 
Oz hues are converted to Ma and Mb structures as 
shown by the dashed (red) lines. There are no true phase 
transitions when crossing these non-solid lines. 

The phase behavior in the inner (grayed-out region) 
of Fig. [njb) is sufhciently complicated that we chosen to 
provide a separate Fig. d] to describe the behavior there. 
The six panels of Fig. [3] show a blow-up of the phase 
diagram in the range Dz £ [0, 0.05] a.u., with color coding 
as explained in the caption. We shall not describe all the 
details here, as these delicate transitions occur in quite 
a small region around the origin in D space and are not 
very relevant to the broader discussion. 

In Figs. [2Ia-d), the solid outer boundary (second or- 
der transition to the rotationless phase) expands to larger 
Dx and Dy with increasing Dz (note the change of scale 
from the first three to the last three panels). In Fig.[3Jc), 
which is for Dz=0.08a.u., the first-order boundaries, 
shown by solid lines marked by squares (magenta), are 
the remnant of the first-order boundaries of Figs. [H^e-f); 
these diminish and disappear as Dz is increased further. 
Then, by the time -D^=0.16a.u. is reached in Fig.[3l^d), a 
new pocket of rotationless phase appears near the origin 
in the D^-Dy plane. This pocket grows until, at a critical 
value of 152=0.172 a.u. shown in Fig.[3]Je), the inner and 
outer regions connect and split the region of rotational 
phases into four ellipses, as shown for Dz—0.2Aa.vL. shown 
in Fig. [Sff). We expected these ellipses to shrink and 
disappear with a further increase of Dz, but in fact this 
happens only very slowly; along the line D^ = Dy = Dz , 
the rotations survive to quite large values of D, as is con- 
firmed by the first-principles calculations upon which the 
model is based. We comment on this further is Sec. lIIIEl 



2. Phase diagram for PbTiOs 

The situation is simpler for the rotational phase dia- 
gram of PbTiOa. The phase diagrams for several snap- 
shots at increasing Dz are shown in Fig. [5] for PbTiOs 
using the same conventions as in Fig. [3]wherever possible. 
When Dz=0 as in Fig.[nia), the center point (D=0) is in 
the TZ phase. The surrounding area enclosed by the solid 
lines (green) is triclinic except along the [100] and [110] 
symmetry lines, which are Max (hues with circles, red) 
and Mbz (dashed lines, red) respectively. The transition 
is continuous across the latter, but first-order across the 
former. Essentially, at small three-dimensional D, the 
system prefers to be in a slightly perturbed version of 
one of four TZ phases, depending on the octant in which 

D resides. Using the notation denotes the octant 

with Dx > 0, Dy > 0, Dz < 0, etc., we find that the TZ- 
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FIG. 5: Phase diagram for rotational phases of PbTiOs under applied D field. Dx-Dy are cuts plotted at (a) 73^=0.00 a.u.; 
(b) 0.02 a.u.; (c) 0.08 a.u. Conventions are similar to those of Fig. [3] (see text for details). 



like phase with 6 ~ [111] is preferred in iS-f-i— |- and S ; 

9 ~ [ill] is preferred in S y+ and ; 6 ~ [111] is 

preferred in \- and S j ; and ~ [HI] is preferred 

in and S The planes Dx=Q, Dy=0, and Dz—Q 

thus form first-order boundaries in this small-D region, 
appearing as solid lines (labeled with circles, red) in the 
2D plots. It follows that the octahedral rotation can 
be "switched" between these 7<^-like (actually, triclinic) 
phases by a small change of external electric displace- 
ment field. 

The area between the inner (green) and outer (blue) 
sohd lines in Fig. [SJa) is the Mcz phase, which be- 
comes Oz along [110] directions (short-dashed-dotted 
lines, brown) and Tx along [100] direction (dotted lines, 
black). The inner solid (green) lines thus represent 
second-order phase boundaries at which as one 

passes to the outside. 

As was the case for SrTiOs, the high-symmetry phases 
for PTO in Fig. [Sfa) become lower-symmetry phases as 
Dz increases, Figs.[SIb-c). In fact, as soon as Dz > 0, the 
entire region inside the solid (blue) boundary is generi- 
cally triclinic. Special cases occur along the dashed lines 
(red), where the symmetry is A^a or Mb, and along 
the dashed-dotted lines (green), which is A4c- When 
Dz is small enough, as in Fig.ISjb), the first-order phase 
boundaries mentioned above are still visible as the solid 
lines with squares (magenta) near the origin, correspond- 
ing to the Oz phase, but with increasing Dz these shrink 
and then vanish, as shown in Fig. [Sjc) . The solid (blue) 
boundary, outside which the rotational phases disappear, 
can also be seen to shrink with increasing Dz, at first 
slowly and then more rapidly, and to disappear by the 
time Dz reaches 0.10 a.u. 



E. Discussion 

There are quite significant differences between the ro- 
tational phase diagrams for SrTiOa and PbTiOs, as 
shown in Figs.|3]and [5] At small D, the major differences 



arise from the fact that the D=0 ground states are differ- 
ent, namely T and TZ respectively. Thus, small applied D 
fields essentially switch the system between T-like phases 
in SrTiOs, or between 7^- like phases in PbTiOa. 

As D gets larger, the behavior becomes rather com- 
plex, but we can identify an important difference that 
can be traced back to the parameter values of the model. 
Namely, we notice a much more isotropic behavior of 
the outer boundary at which the rotations disappear in 
PbTiOa compared to SrTiOa. In PbTiOa, for example, 
we find that the critical magnitude of D at which the ro- 
tations disappear is ~0.09a.u. and ^0.07 a.u. in the [100] 
and [111] D-space directions respectively. For SrTiOa, 
on the other hand, the corresponding values are ^0.14 
and '^0.40 a.u. respectively. In addition to being larger 
(reflecting the stronger tendency to rotational instability 
in SrTiOa), the anisotropy between [100] and [111] di- 
rections is very much greater, with rotations extending 
much further in D space along the [111] direction. This 
can be understood from the coefficients reported in Table 
I. Restricting ourselves to the case of D = {Dq, Dq, Dq) 
and 6 = (6*0, 6*0, 6*0), the critical displacement Dc can be 
obtained from Eq. [5] as D^ = — /i/ (t -I- 2A -I- k) . From this 
we obtain £1^=0.34 and 0.07 a.u. for SrTiOa and PbTiOs 
respectively. The greatly enhanced anisotropy and large 
value of Dc along this [111] direction can thus be traced 
to the small value of the denominator (r + 2A + k) for 
SrTiOa. 

It is useful to put the above results in perspective re- 
garding the dielectric behavior. SrTiOs is known exper- 
imentally to remain paraelectric down to OK, so that 
the entire space of D-fields should be accessible by vary- 
ing an applied £ field, even if the dielectric constant 
is very large. PbTiOs, on the other hand, is strongly 
ferroelectric, so that the region of small D corresponds 
physically to the saddle point of the multi-well energy 
landscape. This internal-energy landscaped as a func- 
tion of D (without octahedral rotations) was mapped 
out in our previous work2^, where the spontaneously po- 
larized tetragonal ground state occurs at l-D[ooi]l = 0.17 
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a.u. Similarly, the spontaneously polarized states with 
constrained orthorhombic and rhombohedral symmetry 
occur at jD^no] |=0.15 a.u. and |i)[iii] |=0.14 a.u. respec- 
tively. So, we can roughly think of this as a three- 
dimensional "Mexican hat" potential with a radius of 
^0.15 a.u. In comparison, the results presented above 
show that the octahedra rotations disappear for \Dz\ > 
0.09 a.u. and |-D[iio] > 0.11 a.u. Thus, the entire re- 
gion of the interesting rotational phase diagram shown 
in Fig. [5] lies inside the Mexican-hat radius, in the re- 
gion where the crystal is unstable under fixed S (but not 
under fixed D) electric boundary conditions. 

IV. SUMMARY AND CONCLUSION 

In summary, we have investigated the phase tran- 
sitions associated with oxygen octahedral rotations in 
SrTiOs and PbTiOa as a function of a three-dimensional 
applied electric displacement field, first directly from 
first-principles calculations and then also using a fitted 
Landau-Devonshire model. For SrTiOa, the D=0 ground 
state is tetragonal, with degenerate states corresponding 
to the rotation angle lying along one of the three Carte- 
sian axes, and for small D-vectors the ground state is a 
weakly perturbed version of one of these states. Similarly, 
for PbTiOa, the D=0 ground state is rhombohedral, with 
four degenerate states having rotation axis in one of the 
[111] or related directions, and again a small D selects 
and weakly perturbs one of these states. However, as 



the strength of D is increased, we find a quite compli- 
cated phase diagram for each material, with both first- 
and second-order phase boundaries appearing in different 
parts of the diagram. The structure is especially rich for 
the case of SrTiOs. For both materials, the general state 
associated with generic ^ Dy is triclinic, 

but states with higher symmetry tend to arise when D 
itself has higher symmetry. In both materials, the rota- 
tions eventually disappear at sufficiently large values of 
applied D. 

Our work represents one of the first attempts to carry 
out a systematic three-dimensional characterization of 
the interplay between polar and octahedral-rotation de- 
grees of freedom in perovskites of this class. While there 
is no external field that couples directly to the rotational 
degrees of freedom, so that is very difficult to find ways of 
controlling the rotations directly, the present work sug- 
gests that such control should be possible indirectly via 
the application of appropriate electric fields. In any case, 
the observed richness of behavior suggests that there may 
be much more to learn in other materials of this class and 
in more distantly related materials. 
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